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ABSTRACT.  In  (1],  1971,  a transformatiSHi^ 
to  as  the  G-Transformation  was  introduced 
Gray,  T.  A.  At^hi^on  and  G.  V.  McWilliams.  "The 
(r-transformation -was^sTiown  to  be  of  value  in  evaluat- 
ing inproper  integrals  and  was  especially  effective 
in  those  cases  where  the  integrand  behaved  like  a 
linear  combination  of  exponentials  with  real  or  com- 
plex arguments.  These  ideas  have  more  recently 
■Y-(Wi21  1970*77)  been  utilized  to  define  a spectral 

estimator  referred  to  as  the  G-Spectral  Estimator. 
In  this  paper  the  nature  of  the  G-transform  emd  G- 
Spectral  estimator  are  discussed  and  it  is  pointed 
out  that  such  an  estimator  is  most  appropriate  for 
processes  whose  autocorrelations  behave  as  a linear 
combination  of  exponentials  with  real  or  complex 
arguments . 


Definition  2.  We  define 


/ f(x)dx;h,n 


J If  (x)-t-f  (-x)  ]dx!h,n 


From  Definition  2 and  our  previous  remarks  it 
follows  that  if  f(x)-tf(-x)  satisfies  a differential 
such  as  (2)  for  some  n,  then 


/ f(x)dx;h,n 


/ f(x)dx 


Since  our  interest  here  is  spectral  estimation 
and  since  we  will  restrict  ourselves  to  real  pro- 
cesses we  will  henceforth  only  consider  symmetric 
integrands.  All  of  our  results  will  however,  with 
obvious  modifications,  follow  from  Definition  2 
when  this  is  not  the  case,  see  (31 . That  is  the  G- 
transform  method  for  estimating  (or  approximating) 
Fourier  transforms  is  equally  applicedsle  to  nonsym- 
metric  integrands.  We  are  now  in  a position  to  de- 
fine the  G-Spectral_Estimator. 

Definition  3.  Let  R(T;T)  be  an  estimate,  based  on 
a record  of  length  T,  for  the  autocorrelation, 

R(T),  of  real  stochastic  process  with  spectral  den- 
sity S (w) . Then  we  define  the  G-Spectral  estimator 
S(w;T)  by 


Let  us  begin  with  a few  necessary  definitions 
Definition  1.  Let  a be  constant  and  let 


where  F(<«)  is  assumed  finite. 

Then  we  define  the  G-transform  by 

IF(t)  F(t-fh)  •••  1 
f(t)  f(t-th)  •••  ; 


F(t-fnh) 
f (t-tnh) 


G(F(t)ih,nl 


f(t)  f(t-t-h) 


f (t-tnh) 


In  (5)  the  quantities  n,  T^,  h are  constants  which 
must  be  determined  in  order  to  calculate  S(w;T). 

A discussion  of  that  problem  is  beyond  the  scope  of 
this  paper.  A complete  solution  to  the  problem  is 
however  given  in  [3],  where  a computer  program  for 
calculating  S(w;T)  is  given.  In  most  instances 
1 < Tg  <<  T,  h * 1,  and  n is  not  actually  needed. 
Again  see  (3] . The  spectral  estimator  defined  by 
(5)  enjoys  many  properties  which  most  other  methods 
do  not.  We  list  a few,  which  subject  to  mild  con- 
ditions can  be  shown  when  n is  properly  chosen,  and 
R(t)  satisfies  an  equation  such  as  (2). 

i)  lim  S(w;T)  = S(w),* 

ii)  if  R(t;T)  - R(t)  when  0 < T < T < T, 

S(w)  i S(w;T) 

iii)  for  most  cases  T.  + (2n-l)h  « T. 


|f(t-Kn-l)h)  •••  f (t-t(2n-l)h)| 

Properly  defining  G in  the  indeterminate  case,  see 
[11 , it  can  be  shown  that  if  f satisfies  the  differ- 
ential equation 


For  brevity  we  will  limit  ourselves  to  inte- 
grals whose  limits  become  infinite  in  only  the  posi- 
tive direction.  Since  the  G-transformatlon  is  non- 
linear it  is  not  clear  from  (1)  how  it  should  be 
defined  for  integrals  whose  limits  become  infinite 
in  both  directions.  Thus  we  include  the  following 
definition. 


* Here  the  limit  can  be  Interpreted  in  whatever 
sense  is  appropriate. 
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In  essence  property (i)  states  that  when  R(t;T) 
is  consistent  for  R(t),  S(w:T)  is  consistent  for 
S(w).  Property  (ii)  states  that  if  R(t)  satisfies 
an  equation  such  as  (2)  (the  coefficients  of  that 
equation  are  of  course  not  known  nor  are  they  as- 
sumed known)  then  there  is  no  error  in  the  approxi- 
mation due  to  leakage,  i.e.  the  finite  record.  Fi- 
nally property  (iii)  is  listed  simply  to  point  out 
that  the  sample  autocorrelation  will  normally  need 
be  calculated  at  considerable  fewer  lags  in  the  cal- 
culation of  S(w;T)  than  other  methods  which  make  use 
of  R(t!T) . 

We  now  give  three  examples  and  some  comparisons. 
In  the  final  two  examples  we  use  discrete  data  and 
the  trapezoid  approximation  to  the  integrals  in  (1). 
when  this  is  done  the  proper  G-transform  will  yield 
the  discrete  Fourier  transform  of  the  corresponding 
discrete  autocorrelation,  i.e.  the  leakage  error 
will  still  be  eliminated.  A more  complete  discus- 
sion will  be  given  in  [3]. 


Exan^le  1.  Consider  the  continuous  autoregressive 
process  defined  by 

^ + ax  = z(t)  , -»  < t < », 


where  z(t)  is  white  noise. 

In  this  case  RCt)  =»  e “1^1  and  for  t > 0 and  h > 0 


2Ttwxdx 


:h 


2a 

2 2 
a +(2irw) 


= S (w) . 


Example  1 is  of  course  trivial  and  only  meant  to 
demonstrate  the  exactness  property  of  the  G-trans- 
form.  A result  which  is  less  trivial  and  suggested 
by  this  ex2unple  and  our  comments  above  is  the 
f ol lowing . 


Example  2.  Let  z(t),  t e (0+1,...)  be  white  noise 
emd  let  {x(t)}  be  the  discrete  autoregressive  pro- 
cess 


x(t)  - Bx(t-l)  » z(t)  . 

Then 

R(k)  - 81*^1 


and 


S(w) 


2 

1+8  -28  cos  2iiw 


Thirty  realizations  of  100  observations  were  gener- 
ated from  this  process  and  S (w)  was  estimated  by  the 
Bartlett,  Tukey  and  Parzen  Window  methods*.  Burg's 
MEM,  Parzen ' s Autoregressive  Spectral  estimator 
(ARSPEC)  and  the  G-Spectral  estimator.  The  summary 
statistics  are  given  below.  The  numbers  in  paren- 
theses show  the  range  of  the  number  of  lags  for 
which  the  seunple  autocorrelation  was  needed  except 
in  Burg  method  where  it  is  the  range  of  the  lengths 
of  the  prediction  error  filter. 


* The  window  methods  here  and  below  are  better 
than  would  be  expected  in  practice  since  several  dif- 
ferent length  windows  were  used  and  the  best  result 
is  given  here. 


SUMMARY  STATISTICS 


(w) 

S(w) 

Estimator 

Sample 

Mean 

Sample 

MSE 

0.00 

19.000 

♦BARTLETT (15) 

7.731 

129.560 

TUKEY (13) 

7.724 

128.775 

PARZEN (19) 

7.964 

123.953 

ARSPEC (1-8) 

14.514 

85.275 

MEM(l-]0) 

15.152 

76.627 

G(2-6) 

15.886 

57.839 

0.06 

1.393 

BARTLETT 

2.684 

2.078 

TUKEY 

3.173 

3.328 

PARZEN 

2.999 

2.838 

ARSPEC 

1.629 

.309 

MEM 

1.764 

.856 

G 

2.146 

2.031 

0.24 

.112 

BARTLETT 

.230 

.0230 

TUKEY 

.162 

.0109 

PARZEN 

.127 

.0109 

ARSPEC 

.131 

.0019 

MEM 

.130 

.0041 

G 

.143 

.0069 

0.50 

.053 

BARTLETT 

.104 

.0042 

TUKEY 

.070 

.0019 

PARZEN 

.099 

.0036 

ARSPEC 

.064 

.0006 

MEM 

.062 

.0012 

G 

.070 

.0014 

Example 

3.  Consider  the  ARMA 

(3,2)  process 

*t 

= l.sx^. 

■1  - ^-=*^^-2  * 

.455Xt,3  + 

"t 

Thirty  realizations  of  length  400  points  were  gen- 
erated from  the  above  process  and  the  various  spec- 
tral estimators  were  again  calculated.  The  summary 
statistics  are  as  follows. 


SUMMARY  STATISTICS 


(w) 

S(w) 

Estimator 

Sample 

Mean 

Sample 

MSE 

0.00 

.017 

BARTLETT  (38) 

.065 

.00238 

TUKEY (33) 

.023 

.00014 

PARZEN(46) 

.018 

.00028 

ARSPEC (7-18) 

.023 

.00015 

MEM (8-20) 

.014 

.00004 

G(4-14) 

.018 

.00019 

0.16 

4.063 

BARTLETT 

3.438 

.6718 

TUKEY 

3.615 

.4314 

PARZEN 

3.598 

.4619 

ARSPEC 

3.907 

1.1851 

MEM 

3.771 

1.2200 

G 

4.022 

.1924 

0.18 

4.480 

BARTLETT 

4.079 

.8192 

TUKEY 

4.226 

.6353 

PARZEN 

4.224 

.6710 

ARSPEC 

5.142 

3.7695 

MEM 

4.624 

2.5005 

G 

4.815 

.7261 

0.50 

.319 

BARTLETT 

.366 

.0172 

TUKEY 

.349 

.0154 

PARZEN 

.347 

.0158 

ARSPEC 

.376 

.0231 

MEM 

.361 

.0252 

G 

.336 

.0083 

1 
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The  above  exan^le  demonstrates  that  the  G- 
Spectral  estimator  tends  to  have  a smaller  mean 
square  error  (M£E)  than  the  other  methods  (this  is 
not  clear  in  the  table  for  the  window  methods  since 
the  best  window  was  used  which  of  course  could  net 
be  done  in  practice) . In  some  cases  it  is  signifi- 
cantly smaller,  in  others  it  is  not.  More  exten- 
sive exair5>les  however  do  demonstrate  that  in  many 
cases  the  G-spectral  estimate  is  significantly  bet- 
ter than  the  others  considered  here  and  hence 
should  be  a us€*ful  addition  to  existing  methods. 

For  further  ex2U(^les  the  reader  is  referred  to  [3] . 


As  a final  comment  it  should  be  remarked  that 
the  application  of  the  G-spectral  estimator  is  not 
limited  to  processes  whose  autocorrelation  satisfy 
an  equation  such  as  (2) . Such  processes  were  only 
considered  here  since  theoretical  results  suggest 
taat  this  is  the  setting  im  which  the  method  should 
be  most  effective.  A more  extensive  Investigation 
would  of  course  be  interesting  and  hopefully  forth- 
coming in  the  future. 
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